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Abstract 

The controversy whether or not the Kardar-Parisi-Zhang (KPZ) equation has 
an upper critical dimension (UCD) is going on for quite a long time. Some 
approximate integral equations for the two-point function served as an indication for 
the existence of a UCD, by obtaining a dimension, above which the equation does not 
have a strong coupling solution. A surprising aspect of these studies, however, is that 
various authors that considered the same equation produced large variations in the 
UCD. This caused some doubts concerning the existence of a UCD. Here we revisit 
these calculations, describe the reason for such large variations in the results of 
identical calculations, show by a large-J asymptotic expansion that indeed there exist 
a UCD and then obtain it numerically by properly defining the integrals involved. 
Since many difficult problems in condensed matter physics of non-linear nature are 
handled with mode-coupling and self-consistent theories, this work might also 
contribute to other researchers working on a large class of different problems that 
might run into the same inconsistencies. 



The KPZ equation [1] for surface growth under ballistic deposition was 
introduced as an extension of the Edwards-Wilkinson theory [2]. The interest in the 
KPZ equation exceeds far beyond the interest in evolving surfaces because of the 
following reasons: (a) The KPZ system is known to be equivalent to a number of very 
different physical systems. Examples are the directed polymer in a random medium 
and the Schrodinger equation (in imaginary time) for a particle in the presence of a 
potential that is random in space and time, (b) The second reason, that is more 
important, to our mind, is that it serves as a relatively simple prototype of non-linear 
stochastic field equations, that are so abundant in condensed matter physics. 

The equation for the height of the surface at the point f and time t, h(r,t) , is 
given by 

fUvV 2 h + ^(Vh) 2 +11 (F,t), (1) 
at 2 

where r|(?,t) is a noise term such that 

(r|(?,t)) = and (r|(f,t)r|(r ',t)) = 2D8(? -f ')o(t-f) , (2) 

so that the constant deposition rate is removed. 

One of the quantities of interest is the roughness exponent a, that characterizes 
the surface in steady state. It is defined by 

([h(r)-h(r')] 2 )oc|r-ff . (3) 

It is well known that for dimensions d < 2 the surface is always rough 
(0 < a < 1), and the two point function <|h q h_ q ^ (h q being the Fourier transform of 

h(r ) ) cannot be obtained perturbatively. Above two dimensions, that is known as the 

lower critical dimension (LCD), the picture changes and there is a transition 
depending on the physical parameters, from a non-perturbative strong coupling 
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regime, that is rough, to a weak coupling regime, characterized by a flat surface, that 
can be obtained perturbatively, and its leading behavior is given by the 
Edwards-Wilkinson model. 

There is however a long lasting controversy regarding the existence of an 
upper critical dimension (UCD) above which the surface is always flat [3]. Various 
approaches, applied either to the KPZ equation directly or to the equivalent directed 
polymer problem, suggest that a UCD does not exist. An important indication for that 
comes from numerical simulations [4] which yields a > at least up to d = 7, well 
above any value obtained for the UCD in the literature. This result is also 
corroborated by two real-space renormalization group approaches [5, 6], that predict 
no UCD as well. On the other hand, finite UCD's are obtained by various field 
theoretical treatments [7-13], and the lid expansion of Cook and Derrida [14]. The 
first to obtain values for the UCD from approaches that gained credibility, by being 
able to give good values for a for d = 2, were Bouchaud and Cates (BC) [11]. They 
considered the equation obtained by Schwartz and Edwards (SE) [8] for the roughness 
exponent and obtained a UCD above which there is no strong-coupling solution. The 
value they give for the UCD arising from the SE method is 2.78. The value they 
obtain for the UCD from their own mode-coupling equation is 2.85 (in an erratum 
published a few months later, they replaced these UCD's with 3.25 and 3.75 
respectively [15]). A more recent result worthwhile mentioning is that of Lassig & 
Kinzelbach [3] that predicted d c = 4 , and in addition were able to prove on a quite 

general assumption that anyway d c < 4 . 

Motivation for the present paper lies in the article of Blum and McKane [10] 
that recalculated the upper critical dimension and obtained using the SE method that 
the UCD is 3.2 and using the BC method that the UCD is 3.6. Being unaware at the 
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time of the erratum of BC [15], we believed that such discrepancies, cannot be 
attributed to slight numerical inaccuracies, and that led us to some numerical 
experimentation with both equations, and we were surprised to discover that indeed 
the variations in UCD's obtained by using identical methods may be large. This 
suggested at once three questions: (a) What is the reason for these discrepancies? (b) 
Can the results of a finite UCD be trusted at all? And If the answer to the last question 
is positive, (c) what are the correct values predicted by those methods? It turns out 
that the above questions and the conclusions that follow are relevant even in view of 
the erratum of BC. In fact this discussion extends well beyond the KPZ equation to 
the growing class of mode-coupling and self-consistent theories, applied to many 
non-linear problems. In addition to answering the above questions we had in mind the 
need to establish rigorous foundations for these theories in order not to face again 
these misfortunate inconsistencies, thus to enable easier development for them. 
Our conclusions to be described later in detail are the following: 

(a) One possible reason for the discrepancies is that both conclusions (BC and 
SE) involve the evaluation of two-dimensional integrals that are not absolutely 
convergent. Therefore, different ways of evaluating the integrals yield different 
results. 

(b) This can be corrected by introducing a high momentum cutoff in the 
original KPZ equation, and letting eventually this cutoff to tend to infinity. This 
defines a unique way of evaluating the integrals that may be described as physically 
correct. 

(c) Even if the correct integration scheme is used, there is still an interesting 
difficulty. The upper critical dimension should be the dimension where the roughness 
exponent, a, becomes zero. It turns out, however, that the integrals under 
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consideration are not continuous as a function of a. If we denote such an integral by 
l(d, a) we find that 

liml(d,cc)*l(d,0). (4) 

This immediately implies that in the equation that should determine the upper critical 
dimension (which is an equation for d) a is not to be set to zero but the limit a — > + 
has to be taken. 

(d) Taking into account the previous remarks we find that indeed the SE and BC 
method yield upper critical dimensions that are 3.7 and 4.3 respectively. This is 
certainly a big difference in the predicted d c 's . In fact, it is so big that while in refs. 

[10] and [15] the BC results obey the Lassig inequality d c < 4 [3], the proper 

evaluation of the integrals in our paper shows that the critical dimension predicted by 
the BC equations violates this inequality. Thus, the proper evaluation of the integrals 
shows that the BC equations are not consistent with Lassig's result, while the SE 
result is. 

In the following we detail the reasons for our conclusions. 
In the SE method [8-9], the roughness exponent a is determined by solving 
the transcendental equation F(d,a) = 0, where 



F(d,a) = -fd d t ? X ^tA J t.e.t-^ + (e-t).e.| e-tt (d+2a) 

J [ t-+| e-t| 2 -%l]L 1 ' ' ' 



+ 



fd d t r [ ''^ % 2 .t- (d+2a) .| e-t| 
J [ t 2 -* + | e-tp + l] ' ' 



+ 

(5) 



(d+2a) 



where e is a unit vector in an arbitrary direction and the t integration is over all 
d-dimensional space. 



denominator 



The BC equation [11] is similarly G(d,a) = 0, where in the definition of G the 

, appearing in the definition of F, is replaced by 



t 2 - a +ie-tr +1 



\2-a 



t 2 ~ a + e-t 



2-a 



From the definition of a it is clear that a must be less than 1. It turns out that 
for such a's the first two integrals in eq. (5) (defined by the two terms in the square 
brackets in the numerator) do not converge absolutely. It is known that an integral, 
which does converge absolutely, can be summed up in different ways, resulting in 
different values. This ambiguity calls for a clear and natural definition of how exactly 
this integral should be performed. 

A physically natural way of defining the integrals is to assume that the original 
KPZ equation has a high momentum cutoff corresponding to a finite size of the grains 
poured onto the surface. The integrals have to be performed under this assumption 
and the size of the grain is then taken to zero. 

The effect of a high momentum cutoff is to multiply the integrand in the 

integrals on the right hand side of eq. (5) by (R - 1) • (R - e - 1 1) , where R is the 

parameter that has to be taken eventually to infinity (clearly a soft rather than an 
abrupt cutoff leads to the same results). It is amusing to consider the two contributions 
to the first integral on the right hand side of eq. (5). These should be identical, 
because one is transformed into the other by the transformation t <-> e - t . The actual 
values of the two contributions depend on the way the integrals are performed. If we 
introduce the cutoff as described above the two contributions are identical, but if we 
do the angular integration first and then do the t-integration from to R, letting 
eventually R tend to infinity (In fact, it can be shown that in calculating I[ or J : we 
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use only 0(R - 1) as a cutoff factor in both integrals. It can be shown that in these 
cases it is enough to use only the factor 0(R-t) as a cutoff), we obtain different 

results for the two contributions. This is not to say that the natural cutoff procedure is 
the only one that yields identical results for those two contributions. 

In order to continue our discussion, let us now apply the above definition for 
the integration procedure, and use a more convenient form for the function F(d,a) of 
the SE equation [8, 9] by defining 



-(d+2a) 



I l( d,a),{^M^I^k 

[t 2 - a + \e-ir + l] 



(6) 



and 



I 2 (d,a)=U V r /J 1 L 1 (7 ) 

[t +\e — t\ +1J 



we arrive at the equivalent SE equation 

F(d,a) = -2I 1 (d,a) + I 2 (d,a) =0. (8) 
Similarly, by defining J ; (d, a) (i = 1,2) the same as I ; (d, a), except that the term in the 



denominator [t +\e-t\ +1J is replaced by [t +\e-t\ J, we can 



reformulate the BC equation [11] as 

G(d,a) = -2J 1 (d,a) + J 2 (d,a) = 0, (9) 
and the four additional schemes for determining a in the article of Blum and McKane 
[10] as 

(scheme 1) d a [-2I,(d,a) + l 2 (d,a)]=0 , (10) 

(scheme 2) d a [-41^,00 + J 2 (d,a)] = 0, (11) 
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(scheme 3) 3 a [I 2 (d,a)] = , (12) 

(scheme 4) 3 a [J 2 (d,a)] = 0. (13) 

where d a means differentiating with respect to a. 

Suppose next that the SE equation does not yield an upper critical dimension. 
In such a case we would expect a roughness exponent a that vanishes as the 
dimension d tends to infinity. (This expectation is consistent with numerical 
simulations as well as with almost all the approaches predicting no upper critical 
dimension, that usually predict oc~ 1/d.) Consequently, the SE equation (eq. 8) for 
large d and small a must take the form 

f 1 (d,a)-f 2 (d) = 0, (14) 

f 2 (d) , , 

where for a> lim — - = and f, (d,a) vanishes when a is set to be zero. 

o-f^a) lV } 

Therefore, we derive in appendix A the large d asymptotic expansion for the integrals 
Ij (d, a) , I 2 (d, a) , J! (d, a) and J 2 (d, a) . 

In appendix A we also obtain the leading order dependence for I, (d,0) and I 2 (d,0) . 



We find that f, (d,0) = -21, (d,0) + 1 2 (d,0) = _ a (^ ^ . It is thus clear that 

the condition for having a small a solution for large d is not obeyed. The case of the 
BC equation is a bit more interesting, because in this case we find that the leading 
order contribution -2J, (d,oc)+ J 2 (d,a) vanishes when a = 0. In order to see whether 

a small a solution exists for large d, a more general examination is needed. The 
argument following eq. (14) assumes that f, (d,oc) is continuous in the limit a— > 0+ . 

In fact, if the function (d,a) is not continuous at a = what has to be obeyed is not 



f x (d,0) = but rather lim fj (d,a) = . As shown in appendix B 
lirn Jj (d,a) = J[ (d,0) . Therefore, the BC equation does not have a small a solution 
for large d. Indeed, also lim I, (d,a) ^ ^ (d,0) but still lim f l (d,a) does not vanish 

ct — ^ 0+ ct — ^ 0+ 

for the SE case. 

Next, we tried to obtain a solution for a that is no necessarily small. This 

implies equating f, (d,oc) to zero for the SE and BC cases. Since f, (d,oc) is obtained 

explicitly this is straightforward and no solution was found. 

The discussion above suggests two reasons for the appearance of severe 
numerical discrepancies. Taking into account the nature of the numerical difficulties 
we have recalculated the upper critical dimension and found from the SE equation 
d c = 3.7 and from the BC equation d c = 4.3 (a result which does not obey the Lassig 

inequality d c < 4). It is interesting to note that the additional schemes by Blum and 
McKane [10] involve integrals that are more benign than Ij or Jj , that are the cause 
of trouble in the SE and BC cases respectively. Therefore, it is quite natural that we 
recover their corresponding result with no deviations, for the three schemes in which 
they find an upper critical dimension (schemes 2, 3, 4 - eqs. (11)-(13) above). 

Another important motivation for this work was the growing interest in 
mode-coupling and self-consistent-expansion approaches in many fields of non-linear 
science. The success of these approaches to handle problems that are not otherwise 
manageable is of course the important reason for their popularity. In this paper we 
intended to make a contribution for the development of these approaches by laying 
more rigorous foundations for them. We showed that in the KPZ problem the lack of 
such a rigorous understanding led to the publication of contradicting results for the 
values of the upper critical dimension as well as for the critical exponents (in high 
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dimensions) by researchers using exactly the same equations. Such a situation 
indicates the need of a better understanding of the equations derived from these 
approaches, an understanding we offer. We believe that once this issue is settled a 
more consistent, fruitful and rapid development of the mode-coupling and 
self-consistent approaches will be possible. 

We thank Sam Edwards for interesting discussions, and Mike Cates for his 
helpful comments, and for drawing our attention to the existence of the erratum [15] 
that yields a higher UCD. 
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Appendix A 

In this appendix we derive the large d asymptotic expansion for the integrals I, (d,oc), 
I 2 (d,a), J[(d,a) and J 2 (d,a). We begin with I,(d,a) and I 2 (d,a). First, it can 
easily shown that the d-dimensional integration in I, (d,a) (eq. (6)) can be reduced to 
the following 2-dimensional integral (using hyperspherical coordinates) 



. rt-(e-t)lrt-elr (d+2<l) 
I (d.oc) = Jd ! t 1 1 ' 



t 2 ~ a +|e-t| +1 

d2 rtcose-t 2 l[tco S e]r (d+2a) 
= n(d-i) t d l dtf(sine) d 2 de^- JL ! 



il-a/2 



(Al) 



■ r> „ il-a/2 

o o i l ~ a + t 2 -2tcos9 + l +1 



Where Q.(d - 1) is the surface area of the unit sphere in (d-l)-dimensions. 
After substituting x = cos 9 we get 



I, (d,cc) = Q(d - l)Jdt Jdx(l - x 2 f 



(x-tjxt 1 " 



2a 



t 2-a + 



t -2tx + l 



l-a/2 



(A2) 



+ 1 



And similarly for I 2 



i d 
I 2 (d,a) = Q(d-l)JdtJdx (l-x 2 ) 

-1 



(x-t) 2 t 12a (t 2 -2tx + l) 



(d/2+a) 



.2-a 



+ t -2tx + 1 



il-a/2 



(A3) 



+ 1 



To evaluate the integrals for large d we expand the d-dependent parts of the 
integrands and perform the x- integration (Laplace method). This gives to a leading 
order in d 
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,l-2a 



I l(d ,a) = a(d-l)^Jdt— - f , 

d o t 2 " + (t 2 + l) +1 



l-(2-a)^ 



(t 2 + l) 



l-a/2 



+ 1 t 2 « + (t 2 +1) 1 01/2 + 1 



(A4) 



^(d-l)^ Kl (a) 
d 



To evaluate I 2 , we note that to apply the Laplace method we have to find the 
maximum of the product (l-x 2 )(t 2 -2tx + l) as a function of x in the range 
-1 < x < 1 and for fixed t. It turns out that for < t < 1 the maximum is attained at 

x = t and for 1 < t < <*> the maximum is attained at x = 1/t . Therefore, we split Jdt 



1 CO 

into two parts Jdt and jdt , where for each of these t regions the x integral is 



performed separately and I 2 is given by 



t " (l — t ) 

I 2 (d,a) = Q(d-l)^ r fdt V ' 

d o t 2 " + (l - 1 2 ) +1 



+ 



+a( d -i)^ dt ' — JL — 

V / jl/2 I t „ / -> \l-a/2 

d i t 2 " + (l - t 2 ) +1 



(A5) 



Notice that in the second contribution to I 2 ,the integral over t still has a d 
dependence. So next we calculate this contribution for large d 



t 1 2a (l-t 2 ) " 
I 2 (d, a) = a (d - 1) ™ jdt 1 > 

d o t 2a +(l-t 2 ) +1 



+ 



+a(d-l)^2 a B(2-a,d-(2-a)) 



(A6) 
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where B (x, y ) is the beta function. 

Now, B(2-a,d-(2-a)) = r(2-a)— j— — — , therefore for large d we use 

v / 

the Stirling formula to obtain the following d dependence of I 2 



I 2 (d,a) = Q(d-l)pjdt - { > 

d o t 2 " + (l - 1 2 ) +1 



l-a/2 * 



(A7) 



4fl(d-l)^2-r(2-a) 



Since a < 1 , the leading order of I 2 (d,oc) is given by 



I 2 (d,a) = n(d-l)^K 2 (a), 
d 



where 



K '(«) = J& 2n / 2 v-^ • (A8) 

o t 2 - a +(l-t 2 ) +1 



The coefficients Kj (a) and K 2 (a) can be easily evaluated at a = to yield 



I 1 (d,a) = a(d-1)^ (A9) 



I 2 (d,a) = Q(d-l)^. (A10) 



Applying exactly the same procedure to J[ and J 2 yields to leading order in 1/d, 
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J l (d,a) = fl(d-l)^L l (a), (All) 
d 



where L ; (a) is obtained from K ; (a) by replacing the expression 



t 2 « + ( t 2 + 1) +1 



in (A4) by 



t 2 « + (t 2 + 1) 



l-a/2 



and 



J 2 (d,a) = n(d-l)^L 2 (a), (A12) 
d 



where again L 2 is obtained from K 2 by the same replacement. 

The coefficients Lj (a) and L 2 (a) can now be evaluated at a = to yield 



J 1 (d,a) = a(d-1)^, 



(A13) 



and 



J 2 (d,cc) = C2(d-l) 



2d 



3/2 



(A14) 
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Appendix B 

In this appendix we show that 



limJ 1 (d,cc)*J 1 (d,0). (Bl) 



a->0+ 



We begin by considering Lj (a) defined in (All) 



L, (a)=}dt — 

J . 2-a 



o t 2 - a +(t 2 +l) 



l-a/2 



1- 



v 2 y 



2(t 2 +iy 



-all 



t +l t 2 - a + (t 2 +l) 



l-a/2 



(B2) 



and we define j(t,a) to be the integrand in this integral. Now, we want to calculate 



limL,(a)= lim fj(t,cc)dt. (B3) 

o-)0+ v ' a->0+ J v ' 



Problems in evaluating the limit can arise from the tail of the integral. Therefore, we 
break the integral into two parts 

R 

Jj(t,oc)dt = Jj(t,a)dt + Jj(t,a)dt (B4) 

R 



where R can be arbitrarily large, and calculate each part separately. 

R 

the first part, i.e. j(t,a)dt, we can take the lim into the integral because the 

J ' <x->0+ 



interval, over which the integration is performed, is finite. 
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R R 



MmJj(t,o)dt = Jj(t,0)dt = J^ T 

Z ' 1 T 1 



2(t 2 +l) 



2t 2 + l 2t 2 +l 



, 1 2R 

dt = -. (B5) 

4 2R 2 + 1 



In the second contribution to the integral in eq. (B4) we expand the integrand for large 
t ' s , assuming we chose a large enough R. We obtain 



j(t,oc) = - 



-l-a 


2t 2-a 


2 t 2-a + J t 2 +1 


j -l-a 


a 1 




— + — r 


2 


2 4t 2 



.l-a/2 



2(t 2 +l) 



l-a/2 



t 2 +1 t 2 « + ( t 2 +l) 
( 1 

V 1 J J 



l-a/2 



(B6) 



Performing the t-integration we get 



jj(t,a)dt = -^R a + 



U-a , (2«-l)(«-2)„ 2a 



^ 1 A 



4 8(a+2) 



R + O 



V R 4+ ", 



(B7) 



Now, we can take easily take the lim 

a->0+ 



f 11 

j(t,a)dt — >- + — - + 

l Jy ' °^ 0+ 4 8R 2 



f 1 A 



(B8) 



At this point we remember that R can be taken to be arbitrarily large, so we obtain 



lim f j(t,oc)dt = 

a->0+J 



1 2R 2 ]_ _L 

4 2R 2 +1 4 8R 2 R ^°° 4 4 2 



1 1 1 

->- + - = ■ 



(B9) 
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and therefore lim L, (a) ? Lj (0) . More specifically 



lim L, (a) = 2L, (0) . (BIO) 

o-)0+ v ' v ' 
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